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Abstract 



o 

(N 

^ ' A massive gravity theory was proposed by Visser in the late nineties. 

, This theory, based on a backgroung metric ba/3 and on an usual dynamical 

• metric Qap has the advantage of being free of ghosts as well as disconti- 

' nuities present in other massive theories proposed in the past. In the 

present investigation, the equations of Visser's theory are revisited with 
^ ' ' a particular care on the related conservation laws. It will be shown that 

a multiplicative factor is missing in the graviton tensor originally derived 
I . by Visser, which has no incidence on the weak field approach but becomes 

^jq' important in the strong field regime when, for instance, cosmological ap- 

plications are considered. In this case, contrary to some previous claims 
found in the literature, we conclude that a non-static background metric is 
required in order to obtain a solution able to mimic the ACDM cosmology. 

(N 
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I. INTRODUCTION 

Studies of the luminosity distance as a function of the redshift of type la 
supernovae suggest that the expansion of the universe is presently accelerated 
[1,2]. These data as well as those derived from the space probe WMAP on the 
cosmic microwave background [3] can be explained by the inclusion of the so- 
called "cosmological constant" term in Einstein's equations. Arguments against 
^ I this possibility have been raised in the literature, in particular the "coincidence" 

problem and the interpretation of such a term as the vacuum energy density. 
However, some authors believe that these objections do not represent real dif- 
ficulties for the theory [4]. They claim that the "coincidence" problem is ill 
defined and that the identification of the cosmological constant with the vac- 
uum energy density is probably a mistake! 

Independently of the reality or not of these difficulties for the ACDM cos- 
mology, alternative models have been proposed in the literature. In a first class 
of models, the usual thermodynamic properties of the constituents of the uni- 
verse are modified and the acceleration of the expansion is driven by a negative 
pressure term associated either to particle production (see, for instance, [5] and 
references therein) or to a bulk viscosity term ([6] and references therein). In 
other class of models, the existence of new fields in nature responsible for the 



1 



acceleration are postulated (scalar fields in GR) as well as modifications in the 
Einstein-Hilbert action {f{R) theories) [7-9] or scalar-tensor theories [9-17]. 

In a particular class of theories (massive gravity), the graviton has a small 
but non-zero mass. These theories have a long history and present several 
difhculties. Fierz and Pauli (FP) noticed that the mass term must be quadratic 
for a Lorentz invariant massive spin-2 theory, otherwise a "ghost" appears in the 
spectrum [18]. A major difficulty with the FP theory is the van Dam-Veltman- 
Zakharov discontinuity [19,20]. In other words, in the linearized FP theory the 
extra scalar mode of the graviton does not disappear and remains coupled to 
matter even in the limit of a vanishing graviton mass and, conseqiiently, in 
such a limit the linear equations of general relativity (GR) are not recuperated. 
Moreover, the prediction of the quadratic FP theory for the light bending effect 
differs from GR, practically ruling out the FP theory. However, Vainshtein 
[21] from the analysis of static spherically symmetric solutions, argued that the 
linear FP theory is only valid for distances larger than a certain scale dubbed 
the Vainshtein radius Ry, which goes to infinity as the graviton mass goes to 
zero. For distances less than Ry, around a static spherically symmetric source 
of mass M, the full non-linear strongly coupled theory has to be considered in 
order to recover GR. In the past years, the problem of the continuous matching 
between solutions inside and outside the Vainshtein radius have been extensively 
debated in the literature [22-25]. 

A different approach was proposed by Visser [26], who introduced a back- 
ground metric not subjected to any dynamical equation. The mass term in this 
theory depends both on the dynamical and on the background metric in such 
a way that in the linear limit, the massive field obeys a Klein-Gordon equation 
with a source term. GR is recovered when the graviton mass vanishes. The 
massive field in Visser's theory has six degrees of freedom: five spin-2 and one 
scalar [26]. This theory was applied to cosmology by different authors [27-30], 
who claim that the resulting dynamical equations based on such a theory are 
able to explain the present observed acceleration of the universe and to satisfy 
other cosmological tests like the distance scale provided by the baryon acous- 
tic peak and the cosmic microwave background shift parameter. In the afore- 
mentioned investigations, a flat background (Minkowski space-time) had been 
adopted despite the fact that according to Visser, "in a cosmological setting it 
is no longer obvious that we should use the Minkowski metric as a background". 
If a flat background is adopted, then the only issue resulting in an expanding 
dynamic metric is to assume that the divergence of the massive graviton tensor 
is a source for the divergence of the matter stress-energy tensor [27]. In Visser's 
theory the divergence of the graviton tensor is set to zero and, consequently, 
when the dynamical equations are linearized around the background metric, the 
Hilbert-Lorentz condition appears naturally and not as a gauge one. It is worth 
mentioning that other bimetric gravity theories with generally a flat prior ge- 
ometry have been elaborated in the past, in particular the approach by Rosen 
or the vector-bimetric theory by Rastall (see, for instance, [32,33] for reviews). 
These theories are quite distinct since the lagrangean density from which the 
fleld equations are derived differs drastically from that proposed by Visser. 

In the present paper, the Visser's theory is revisited as well as applications to 
cosmology. We will show that the graviton tensor derived by Visser [26] must be 
corrected by a factor equal to the ratio between the background and dynamical 
metric determinants. The conservation laws are also revisited and the following 
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question is examined: if the divergence of the graviton tensor is equal to zero (as 
it should be expected from general arguments based on a variational approach 
of gravity) is it possible to find an adequate background metric from which 
a cosmology compatible with the present data emerges? A positive answer 
can be given but the solutions able to mimic the ACDM cosmology are not 
completely satisfactory, since they require that the background metric tensor 
be proportional to the scale factor describing the dynamics of the universe. 
This paper is organized as follows: in Section II the Visser's theory is revisited 
and the correct graviton tensor is derived. An alternative definition of the 
lagrangean density describing the graviton field is given as well as the resulting 
field equations. In Section III an application of the theory to cosmology is 
discussed and finally, in Section IV the main conclusions are presented. 

II. THE FIELD EQUATIONS AND THE CONSERVATION LAWS 

In this paper, the notation ruap = diag (—1, +1, +1, +1) (Minkowski metric 
in cartesian coordinates) will be adopted as well as units such G = c — h — 1. 
The action proposed by Visser [26] reads 



S= (Tx 



1 

16^ 



gR + V-bLmass (.9, b) + \/^Lmatter {g,X)> (1) 



where bap represents the background metric tensor, g^p the dynamical metric 
tensor and X stands for any non gravitational field. Both metrics are required to 
have a lorentzian signature (— , +, +, +). The lagrangean of the massive graviton 
is given explicitly by 

L^ass {g, b) = -im2 (^^^ _ (^^^ _ _ 1 (^^^ _ 2 

(2) 

The contravariant tensor 6"^ is defined from bap by the inversion relation 
bapb^'' = ^2: while g"^ is defined from the relation gapg^^ = <52 as usually. 
The coefficient - 1/4 ensures that in the weak field limit, i.e., when b^p — i^ap 
and ga/s = map + hap with \hap\ << 1, the graviton field in vacuum obeys the 
Klein-Gordon equation. Since the construction of the action defined in eq. [T]is 
mainly motivated by such a requirement, it is worth pointing that any alter- 
native action leading to the same linearized equations would be acceptable a 
priori. In particular, the graviton lagrangean density can be defined in terms of 
the dynamical metric, i.e., \J —gLmass instead of sf—bLmass- This corresponds 
to an alternative to Visser's proposal that admits the same usual weak field 
limit. Such an alternative for the graviton lagrangean density will be examined 
in some more detail at the end of this section. 



A. Conservation laws 



From the variation of the action (eq. [1} with respect to the field X (here, 
for simplicity, we assume the presence of only one non-gravitational field) one 
obtains the Lagrange equation describing the dynamics of the considered field. 
Since the matter lagrangean density Lmatter depends only on the field X and on 
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the dynamical metric tensor (the field X couples with gap only), then the diffeo- 
morphism invariance of Lmatter leads immediately to a conservation equation 
expressed by the null covariant divergence condition 

VaT"'' = (3) 

where T""^ is the stress-energy tensor of matter [31]. If the condition above 
is ignored, some inconsistencies may appear in the physical laws describing 
the dynamics of non-gravitational fields. It is worth recalling that the null 
divergence of the Einstein's tensor leads only to the null divergence of the sum 
of tensors constituting the right side of the field equations derived by varying 
the complete action of the theory (in our case eq. [T]) with respect to the metric 
tensor gap- If different stress-energy tensors are present, the null divergence 
of the Einstein tensor is not equivalent to eq. |3l In this case the validity of 
eq. [3] implies an additional condition to be fulfilled by the solutions of the field 
equations. 

B. Field equations 

Varying the action defined by eq. [T] with respect to the metric tensor gap 
leads to 

G^^ =S^T^^+%^T:i,, (4) 

where G"^ = R°'^ — ^Rg"^ is the Einstein tensor. Rap is the Ricci tensor and 
R = R^ is the Ricci scalar. The matter stress-energy tensor T"'^ is, as usually, 
defined by 

j,ap _ 2 5{\/~gLmat) ,r\ 
I ~ Ji^.aP ^ ^ 



-9 

and the so-called graviton tensor by 

,2 



TZLs = -^^ (b'^'-b'^^ - h^-^b-P^ {g,, - h,,) (6) 

From the null divergence of Einstein's tensor and eq. [3l it results 

V.T,r,L, = 0. (7) 

It should be emphasized that both conservation laws expressed by eqs. [3]and[7] 
must be satisfied and these are necessary conditions to be taken into account 
when considering solutions of the field equations. From these equations, it is 
trivial to verify that Visser's theory reduces to GR when the graviton mass m 
vanishes. As we shall see in details in appendix A, a class of solutions in which 
the dynamical metric tensor is proportional to that of the background met- 
ric exists, which can be generated from GR solutions including a cosmological 
constant. 

C. An alternative lagrangean density 

As already mentioned, the requirement that the graviton field obeys the 
Klein-Gordon equation in the weak field (or linear) limit does not fix the form 
of the lagrangean density. An alternative to Visser's proposal is to consider 
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a lagrangean density defined in terms of the dynamical metric tensor, i.e., 
y/^gLmass and, in this case, the total action of the theory is 



S 



y"(i*a;^^|Y^ [R+ Ljnass {9,b)] + Lmatter {g , X)^ 



(8) 



The resulting field equations are essentially identical to eq. |3]but now the gravi- 
ton tensor reads 



771 

16^ 



(9) 



m 
16^ 



It should be emphasized that as in the case of the original Visser's proposal, 
the field equations derived from this alternative lagrangean density admit so- 
lutions generated from GR equations including a cosmological constant (see 
details in appendix A). 

Ill - FRW COSMOLOGY 

Recently, different authors have considered the Visser's theory to describe 
the dynamics of the universe ([27-30]), claiming that the graviton mass term 
appearing in the Friedman equations is able to drive the observed acceleration 
of the expansion of the universe. In these investigations, a Minkowski spacetime 
was assumed as a background and, in this case, as a consequence of eq. [71 it 
results that the scale factor must be constant or, in other words, the universe 
must be "static" . In order to maintain a flat and static background (Minkowski) , 
the aforementioned authors assumed that the conservation law expressed by 
eq. [7] is violated and, in order to satisfy the condition VaG^ = 0, they have 
hypothesized that the divergence of the graviton tensor is a (negative) source 
term for the divergence of the stress-energy tensor of the matter. 

Here the FRW cosmology is revisited in the context of the Visser's theory but 
with preservation of the conservation laws as discussed in the previous Section. 
Considering the standard FRW form for the dynamical metric we have 



ds^ = ~dt^ + a {ty dx^ + Fk (x) 



(10) 



where k = +1,0,-1 corresponds respectively to spherical, flat or hyperbolic 
spatial sections and, accordingly, Ffc (x) = sinx, X: sinhx- Following Visser 
([26]), we will search possible non-static solutions for the background metric. 



da^ 



-B{tfde + A{tf dx^ + Gk{xf {de^ + 



sir? Odip^) 



(11) 



where A and B are positive functions of the cosmic time. From the spatial 
components of eq. [7]we have necessarily Gfc(x) = ^fc(x)- The function A{t) 
may be interpreted as a "background scale factor" while a "background cosmic 
time" may be defined by / B {t) dt. Matter, as usually, is described as a perfect 
fluid with energy density e {t) and pressure P (t) , linked by an equation of state 
P ^P{e). 
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Besides the time component of matter conservation (eq. [3]) , the other equa- 
tions describing the dynamics of the universe are the components (0 0) and (11) 
of the field eq. SI since the remaining components do not provide any additional 
information. On the other side, a relation between the coefficients of the back- 
ground and of the dynamical metric can be obtained from eq. [T] It is useful to 
define the functions E{t) = A{t) /a{t) and {t) = B (t) /S(t), which should 
be positive since A, B and the scale factor a are positive functions. In terms of 
these new functions the background metric can be expressed as 



E {tf {tf dt^ + a {tf \dx^ + Fk ixf (de^ + sin' Odip^)^ } (12) 



The usual procedure in the context of the Visser's theory is to define first a 
background metric and then, from the field equations and conservation laws, to 
derived the dynamics of the universe (supposing that the properties of matter, 
like the equation of state are given a priori) . In this case, any cosmological scale 
parameter a{t) can be obtained by an appropriate choice of the background but 
this not necessarily implies that we have a physically acceptable solution. An- 
other possibility to tackle the problem is to search for metric coefficients of the 
background spacetime once a dynamical solution is imposed to the field equa- 
tions or, in other words, to work in the opposite route of the usual procedure. 
This approach will be adopted here. From the field components (0 0) and (1 1) 
of eq. m we have 

' k 8tt 

Y' ' 12 



2 o^+— S (13) 



and 



a fa 



+4 = -8.p+!=!ei!iii^ ,14, 

a \a J 4 



where the notation X = dX/dt was adopted (It is worth mentioning that in 
above equations S and ^ are functions of the cosmic time t) . The new function 
S in eq. [13] is defined as 

S ^ (1 + 2*2S2 - 3*') (15) 

The space components of eq. [7] are trivially satisfied while the time component 
provides a relation between the coefficients of the background metric and the 
scale factor of the dynamical metric, namely 

f + 3°fE + i±il-^2!2!).0 (16) 

dt a\ * / 



or, equivalently 



dTi — 1 / o dlna 

+2*^)^ (17) 



For the sake of completeness, we mention that the time component of eq. |3]can 
be written explicitly as 

' + 3Pa^d ^ 0. (18) 
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Combining eqs. [T31 E] and [12] one obtains 



^ = -| (e + 3P) + ^ [2*^-^ (3*^ - 1) - (3*4 + 1)] (19) 

The equation above shows that the last term on the right hand side (depending 
on the graviton mass) can give a positive contribution to the acceleration of the 
expansion of the universe. According to our adopted approach, a solution for 
the background metric functions S{t) and ^'(i) can be obtained by adopting the 
following procedure: firstly, from the matter equation of state and eq. 1181 the 
variation of the energy density and of the pressure can be derived as a function of 
the scale parameter a or, equivalently, of the redshift z. Then, if the acceleration 
parameter q(z) = —aajk? and the Hubble parameter -ff(z) are known from 
observations, the metric functions S(z) and 'i'(z) can be determined from eas.[T51 
and [191 Acceptable solutions require that the background metric functions be 
positive and at the present time, in order to have a positive acceleration, the 
following condition must be satisfied (including all the physical constants) 

[2*^S^(3*^ - 1) - (3*4 + 1)] > (20) 

where is the present matter density parameter and the metric functions 
S and * are taken at the present time. The other symbols have their usual 
meaning. 

Although static solutions have only an academic interest, it is worth men- 
tioning that these solutions necessarily implies also a static background geom- 
etry and vice-versa. Hence, the ACDM cosmology cannot be reproduced if a 
Minkowski background is adopted. This point will be considered in some more 
detail in appendix B. 

A. Back to the ACDM cosmology 

The present observations (luminosity distance of type la supernovae, the 
baryon acoustic peak (BAO) and the cosmic microwave background (CMB) shift 
parameter) are quite well fitted by the so-called ACDM model. Is it possible to 
mimic such a cosmological model within the framework of the Visser's theory? 
A positive answer can be obtained if one identifies the last term on the right 
side of eq. [13| as 

In this case the function S must be a constant and eq. [T7| has two possible 
solutions. The first corresponds to ^^=1 and = 1 -|- 2m~^A. Since A > 0, 
we have > 1. This particular solution corresponds to a class of solutions 
discussed in appendix A, referred to the cosmological case. The second possi- 
bility corresponds to = (1 — 25^)~^ and E = — 2Vl — 22^. However, this 
solution implies a negative cosmological constant, which is not supported by 
observations and will be not discussed further. 

B. ACDM and the alternative Visser's theory 
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In the framework of the alternative definition of the graviton lagrangean 
density, the (0 0) component of the field equations is formally identical to ea.[T51 
but now the function T, is defined as 

5 - (18 - 12S2) *2 _ (3 _ 12S2 _^ 8^4) *4 
E = 5^ ^ j^^- ^ 22 

Another equation relating the metric coefficients can be obtained from the con- 
servation law expressed by eq. [7] whose time component is 

dE — 1, „^9„,^dlna 

— =3—— (1-^^ + 2^^E^) ——. 23) 
dt *424 ^ ^ ' dt ^ ^ 

Again, to mimic the ACDM cosmology we require that E=constant and, in this 
case, two solutions are possible: 'i>^=l and = (1 — 2S^)~^ (Notice that the 
solution ^^=1 corresponds to a case discussed in appendix A). Contrary to the 
solutions derived from equations based on the original Visser proposal, here the 
two solutions are compatible with a positive cosmological constant and a lower 
limit for the gravition mass can be derived from the observed value of A. In 
order to have at least one solution with a positive cosmological constant, the 
mass of the graviton must be higher than a critical value given by 

m > 2^/m^^ ~ 7.7 X 10"*'^ g (24) 

where we have adopted 57a = 0.7 and Hq — 70kms~^ Mpc~^ . 



IV - CONCLUSIONS 



In the present investigation the Visser's theory was revisited and, in par- 
ticular, we have found that the original graviton tensor must be corrected by 
a factor equal to the ratio between the square root of the determinant of the 
background metric and that of the dynamical metric. This correction is not 
relevant when the usual weak field approximation is considered but is of fun- 
damental importance in cosmological applications. We have also considered an 
alternative to the graviton lagrangean density proposed by Visser and the con- 
sequent modifications in the graviton tensor but both approaches lead to the 
same linear field equations. 

We have also shown that the field equations of the theory when combined 
to the conservation laws are able to mimic a ACDM cosmology if and only if 
the background metric functions are proportional to the scale factor defining 
the dynamical metric. This is in contradiction with claims in the literature 
based on investigations considering a Minkowski background and an abandon 
of the conservation laws expressed either by eq. Oor eq. [71 However, the results 
are somewhat disappointing since only a particular form of the background 
metric, whose choice was based on a priori cosmological considerations leads to 
satisfactory results. 

Since the background required for compatibility between the theory and 
cosmological observations turns to be time dependent, a natural question raises: 
how this affects locally physical process occurring in the linear regime? As 
argued by Visser, in the usual weak field approximation the background metric 
should be Minkowiskian and, in fact, an expanding background metric can be 
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put in the form niap by a suitable choice of coordinates if the relevant timescales 
are orders of magnitude smaller than the cosmological timescale. This is the 
case when planetary motions or the propagation of light within scales of the 
order of the solar system are considered. However if this is true at a given 
instant, this may not be the case (say) some billions years later. In this situation, 
performing an expansion of the dynamical metric a Klein-Gordon equation is not 
obtained. These considerations may lead into specific cosmological signatures 
on the local dynamics, opening eventual new constraints on the graviton mass 
by local observations. These aspects will be analyzed in a forthcoming paper. 
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APPENDIX A - AGR solutions vs Visser's theory 

For a given background metric 6^^, let us search for solutions in which the 
dynamical metric gi^^, is proportional to the background metric tensor, i.e., 

gfii' = A^^&pi. (25) 

where is a positive constant. In this case, from eqs.[B](see text) one obtains 

r^i. = ^-'ff"'' (26) 
Consequently the field equation takes formally the following structure 

^ ^L^rn^gO^P = 87rr"'5 (27) 

Notice that the above equation satisfies consistently eq. [7]and indicates that all 
GR solutions including a cosmological constant A are also solutions of Visser's 
massive gravity with a background metric tensor 

= (l + il) 9o.p (28) 

This result indicates that solutions based on eq. [25] are possible only if the 
background metric tensor is proportional to a solution of Einstein equations 
including a cosmological constant and with the same matter stress-energy tensor. 

Similarly, if the alternative lagrangean density for the graviton field is con- 
sidered, one obtains for the graviton tensor 

that replaced in the field equation leads to 

Thus, as in Visser's original theory, all GR solutions including a cosmological 
constant term are also solutions of the massive gravity field equations but now 
the proportional constant is given by 

(1 - A-) (2A- - 1) = § (31) 



APPENDIX B : Static cases in Visser's theory 
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Static solutions are those with a constant scale factor a (or a — Q). The 
equations to be considered in this case are eqsHSl [TH and [TOl (see text). From 
these equations one obtains 

^ = (32) 

where 

E = [1 + 2*2^2 - 3*2] (33) 
k 11 

— = -SttP + -m^— \2^^'E? -^^ ~l] (34) 
a2 4 * 

with E = (static case), which implies that E =constant. The constancy of the 
energy density and of the pressure can be derived from the first equation and 
the equation of state. Thus, both E and the quantity [l + vI/2 - 
are constants, indicating that the background metric is necessarily static. 

For a given e, P, k and a (the case P — {) corresponds to a dust filled static 
universe in the context of Visser's theory), the metric coefficient ^E* is obtained 
by solving the equation 

- - Stte U-^ - m^-^^ - 2 — + SttP ) * + = (35) 



and, consequently, S satisfies 



*^S^ = ^-Af4 + 8.P)* (36) 



The equations above permit to obtain the metric coefficients ^ and S. 

Reciprocally, if we have a static background, we have necessarily A and B 
or vPS and Sa constants (Notice that these imply also that a/4' is a constant). 
Inserting these relations into eq. [17] and after a straightforward calculation one 
obtains 

^ ' at 

Thus, one has either 2*2^2 _ \i/2 ^ = o or d^/di = 0. Both possibilities imply 
*=constant and a— constant, i.e., a static universe. 
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